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ABSTRACT 

Isothermal models and other simple smooth models of dark matter halos of gravitational lenses often 
predict a dimensionless time delay Ho At much too small to be comfortable with the observed time delays 
At and the widely accepted Hq value (~ 70 km/s/Mpc). This conflict or crisis of the CDM has been 
highlighted by several recent papers of Kochanek, who claims that the standard value of H favors a 
strangely small halo as compact as the stellar light distribution with an overall nearly Keplerian rotation 
curve. In an earlier paper we argue that this is not necessarily the case, at least in a perfectly symmetrical 
Einstein cross system (Paper I, astro-ph/0209191). Here we introduce a new mass degeneracy of lens 
systems to give a realistic counter example to Kochanek's claims. We fit the time delay and image 
positions in the quadruple image system PCI 115+080. Equally good fits are found between lens models 
with flat vs. Keplerian rotation curves. Time delays in both types of models can be fit with the standard 
value of Hq. We demonstrate that it may still be problematic to constrain the size of lens dark halos 
even if the data image positions are accurately given and the cosmology is precisely specified. 

Subject headings: cosmological parameters — dark matter — distance scale — gravitational lensing 



1. INTRODUCTION 

Gravitational lensing provides a powerful tool to con- 
strain the dark matter halos of galaxies. One of the 
promises of gravitational lenses is to constrain the Hub- 
ble constant. However, this has been hampered to some 
extent by the intrinsic degeneracies in models of the dark 
matter potential of the lens (Williams & Saha 2000; Saha 
2000; Saha & Williams 2001; Zhao & Pronk 2001). Now 
that the value of Hq is fairly well constrained by indepen- 
dent methods, e.g., H = 72 ± 8 km/s/Mpc from the HST 
key project (Freedman et al. 2001), and the cosmological 
model has been determined at more and more precision, 
it is interesting to ask whether we can reverse the game 
and set more stringent constraint on the dark matter po- 
tential. To this end, we would like to understand whether 
the Hubble constant and the lensed images could uniquely 
specify the dark matter content, or whether there are very 
different lens models with identical Hq value. 

It is well-known that isothermal models and other sim- 
ple smooth models of dark matter halos of gravitational 
lenses often predict a dimensionless time delay HoAt much 
too small to be comfortable with the observed time delays 
At and the widely accepted Hq value (~ 70 km/s/Mpc). 
Models with isothermal dark halos tend to yield an Hq 
around 50 km/s/Mpc. This conflict has been highlighted 
by several recent papers (Kochanek 2002a, b,c). Kochanek 
(2002a) found that it is difficult to reconcile the time de- 
lays measured for five simple and well-observed gravita- 
tional lenses with Hq ~ 70 km/s/Mpc unless the lens 
galaxy has a nearly Keplerian rotation curve with the halo 
following the stellar mass profile by a constant mass-to- 
light (M/L) ratio. If the lenses had a more plausible flat 
rotation curve (isothermal mass distributions) he found 
Hq = 48^4 km/s/Mpc, which is grossly inconsistent with 
the HST Key Project. Kochanek (2002c) argued that more 



realistic models with a CDM halo plus adiabatically cooled 
baryons behave like isothermal models. They produced a 
still too low Hq unless one adopts a problematically high 
baryon fraction fi(,/f2 m > 0.2 of the universe, and require 
all these baryons to cool. His conclusion was that either 
Hq ~ 70 km/s/Mpc is too high or any lens mass models 
for the observed time delay systems must follow a compact 
distribution, nearly like that of the stellar light, hence has 
very little extended dark matter halo. This argued for 
the first time a new problem for dark matter halos, and a 
particularly serious problem for current CDM paradigm of 
galaxy formation. 

Here we discuss the effect of a new degeneracy in strong 
lensing models in resolving this new problem. It is shown 
that the observed time delays and image positions cannot 
uniquely determine the extent of the lens mass distribu- 
tion. In particular, a system with a very extended dark 
matter distribution could minic a system without any dark 
matter as far as strong lensing data are concerned. Here 
we give an analytical explicit illustration of the degree of 
degeneracy in lens models. 



2. MODELS 

Consider fitting a general quadruple image system with 
the images at cyclindrical radius (Ri,0i) for i = 1,2,3,4 
from the lens center; we shall call the radius to the out- 
ermost image the Einstein radius Re- These images lie 
at extrema points of the time delay surface. All lensing 
properties can be derived from the time delay surfaces. 
As a minimal model, let's consider a spherical lens po- 
tential (f>(R) = <j>*(R) + <ph(R) for the stars and the halo 
plus an external shear potential <j) e (X,Y), the time de- 
lay contours are determined by a dimensionless time delay 
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t(X, Y) given by 



T = 



100/i zi, z s )day arcsec -2 
(X - X s ) 2 + (Y — Y s f A 



(1) 



(2) 



where R = V 'X 2 + Y 2 is the cylindrical radius, all in units 
of arcsec, and ((ft, zi,z s ) ~ 1 is a constant containing all 
the dependence on the cosmology f2 and lens and source 
rcdshifts zi , z s . 

A minimal model for the external shear is 

KJCY)- *"-^**™ . (3) 

A minimal model for the stars and halo would be a spher- 
ical lensing potential 
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where (1 — i/)mo is the total stellar mass enclosed, R e is 
the half-mass radius, and 2 — a specifies the cuspiness of 
the stellar distribution. Here the halo mass distribution 
follows the mass of the stars with an adjustable parameter 
< v < 1, apart from a to-be-determined halo component 
4* U (R). The latter is given by 



4> U {R) — rriQvN 
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where Ri are respectively the radii of the images i = 
1, 2, 3, 4, and the function N(r) is given by 



N(r) = - In [l+ (r 2 - l)' 
8 - 



if r < 1, 
- 1, 



(7) 

(8) 
(9) 
(10) 



- 0.0866 (1 - 2r 2 ) 

~0[(r-l) 3 ] ifr 

<~ lnr if r 3> 1. 

Figure 1 illustrates the function N 4 (r) and the first and 
second derivatives of N(r). 

The potential <p v is designed to have the following prop- 
erty: It dominates the stellar potential <j>* (R) at radii much 
larger than the Einstein radius Re with 

<t> v ~ m a v In 4 (RI R E ), cj>*(R) - (1 - v)m ln(R/R E ). 

(11) 

But at small radius <fr n u is small and has virtually no effect 
on lensing. In particular, it keeps the convergence k, shear 
7, the components of the amplification matrix \i and the 
time delay between images unchanged at the radii of the 
four images R = Ri. In fact it has a vanishing contribution 
to the arrival time surface near the images up to the order 
0(R — Ri) 3 . This can be verified by Taylor expansion. It 
is easy to show (also cf. Fig. 1) that the factor N(r) and 
its first and second radial derivatives obey 
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(12) 



\R=Ri - ly \R=Ri 

Becaue (j> v {R) is a multiplications of the four N(r), all of 
the zeroth, first and second derivatives vanish at the four 
radii Ri. And because <j> v (R) is azimuthally symmetric, 
all its azimuthal derivatives with respect to 9 also vanish. 



The time delay surface 

(X - X s ) 2 + (Y- Y s ) 2 
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Hence when we vary the parameter v, the time arrival 
surface t(X,Y) yields the same extrema, or images. This 
means models with different v (or halo) will have good- 
ness of fit to all strong lensing data (the time delay, the 
image positions and even the amplifications and parity of 
the images). It only alters the mass distribution of the 
lens, e.g., the total mass of the lens and the spatial extent 
of the lens. Hence it creates a degeneracy in lens modeling, 
making it problematic to draw unique conclusions on lens 
halo mass from image modeling. We note, however, by 
construction this degeneracy applies only to point images 
or point sources. 

The surface density and mass of the model can be com- 
puted as 



k(R) = 1 - -V 2 r, 



M{R) = R 



dR' 



(14) 



in particular, the density and the mass for the stars only 
are given by 

a(l — is)moRe 



k*(R) = 



and 



M*(R) = 



2R 2 - a (R a + R*) 2 
(1 - v)m Q R a 



(15) 



(16) 



R a + R% 

As we can see the mass in stars has a half-mass radius 
R = i? e , and converges to a finite mass mo at infinity. 
The density has an inner cusp 2 — a, and drops steeply 
with radius, by a factor of 25/4 from R — R e to R = 2R e 
for cored models (a = 2) or by a factor of 4 for isothermal 
cusped models (a = 0). 

The halo component 4> U (R) contributes to k(0) by a 
negligible amount ±0.001^ for < R < 1.3R E , but con- 
tributes to k(20Re > R > 2Re) by an amount of the 
order vRe/R- This flattens the rotation curve at large 
radii. A flat rotation curve corresponds to a constant de- 
flection strength, 

U =dR^R> 
which is effectively the rotation curve squared. At large 
radius we have 

M -» mQ + 2,vm u \n 3 (R/R E ). (18) 
In comparsion to the case with only the stars, we have a 
nearly Keplerian rotation curve beyond R = 2R e . 

3. RESULTS 

PG1 115+080 is a quadruple system with a nearly ax- 
isymmetric stellar lens at zi = 0.31, and the quasar source 
is at z s = 1.72. We use the photometric data of Impey 
et al. ((1998)) for the images Al, A2, B, C, and the time 
delay tsc = 25 days between the nearest B image and 
the furthest C image (Schechter et al. 1997). There is an 
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infarcd Einstein ring of radius Re ~ 1.4". The stellar lens 
is well-approximated by a de Vaucouleurs profile with a 
half-light radii of R e = 0.55" = OARe- Since the surface 
density is nearly cored, we set the cusp slope 2 — a = 0. We 
adopt a ACDM cosmology with the time delay constant for 
at the redshifts zi = 0.31, z s = 1.72 of PG1115+080 given 
by 

100/Iq ^(^ = 1 - A) ~ (30 - 32)^ May arcsec" 2 . (19) 

Let's try to fit PG1115 with a halo-free model with 
v = 0. Totally we have five free fitting parameters, with 
(71, 72, mo) for the lens model, and (X s , Y s ) being a pair of 
coordinates for the source. Since the four images provide 
in general at least eight constraints, we do not expect a 
perfect fit. Nonetheless a reasonable fit can be found for 
this star-only model, as shown by the time delay contours 
in Figure 2a and the cuts in Figure 2b. The model yields 
a time delay 

t~bc ~ 0.56 arcsec 2 , (20) 

between images B and C. Compared with measured delay 
tsc = 25 day, we find 

h = „ 3075c „ ?2 (21) 

tBC tBC 

for this star-only model. The high value of ho is related to 
the very small convergence at the Einstein radius Re- 

Now consider adding the component <f> v , i.e., we increase 
the value v from zero to v = 0.04,0.08,0.16. We see no 
detectable differences in the time delay (cf. Fig. 2a and 
Fig. 2b), hence we still have ho — 0.72. All models predict 
amplification ratios among the four images to be 

A x : A 2 : B : C = 3 : 4 : 0.7 : 1, (22) 

independent of v. These are in good agreement with the 
observed infared flux ratios 

A 1 : A 2 : B : C = 3.84 : 2.67 : 0.64 : 1, (23) 

apart from the well-known problematic flux ratio between 
the close pair A\ and A 2 , which has been attributed to ei- 
ther microlensing or lensing by substructures (e.g., Impey 
et al. 1998; Barkana 1997; Metcalf & Zhao 2002). Con- 
sistent with Kochanek (2002b), our model density is small 
near Re with a convergence 

k(R e )~ 0.1 + xv, (24) 

independent of v; by construction </>„ has a vanishing con- 
tribution to the convergence at the images. But by raising 
the parameter v the models develop a very massive halo at 



radii 3> Re with a nearly flat rotation curve (cf. Fig. 3), 
very different from the model with v = 0. This illustrates 
an explicit counter example to Kochanek's claim of a con- 
flict between lens models with an extended dark halo and 
ACDM cosmology with Hq <~ 70. We do not find such a 
conflict. Instead we find that the new mass degeneracy 
prevents us from drawing a robust conclusion about the 
dark halo on the basis of the image positions and time 
delays alone. 

4. SUMMARY 

As we can see, it is possible to construct many very dif- 
ferent models with positive, smooth and monotonic sur- 
face densities to fit the image positions. There are also no 
extra images. These models also fit the same time delay 
and time delay ratios using a Hubble constant and cosmol- 
ogy consistent with ACDM cosmology. Hence the models 
are truely indistinguishable for lensing data. They fit the 
flux ratios equally well, and produce nearly indistinguis- 
able Einstein rings, which is the region of minimal gradient 
of the time delay surface. The models have identical light 
profiles, undistinguishable by data. 

Among the acceptable models to PG1115+080, there are 
models with a Kcplcrian rotation curve and models with a 
nearly flat rotation curve. So lensing data plus Ho cannot 
uniquely specify the mass-to-light ratio of this system. 

We conclude that strong lensing data may not uniquely 
determine the Hubble constant, even if we fix the cosmol- 
ogy, the lens and source redshifts, and the time delays 
and amplification ratios of the four images. There are 
at least important degeneracies in inverting the data of a 
perfect Einstein cross to the lens models and the Hubble 
constant. The relation between the value of Ho and the 
size of the halo is not straightforward: a high Ho does not 
necessarily mean no dark halo, and models with a flat ro- 
tation curve do not always yield a small H . We also com- 
ment that it would be difficult to determine the cosmology 
from strong lensing data alone because the non-uniqueness 
in the lens models implies that the combined parameter 
h^^fl, zi, z s ) is poorly constrained by the lensing data, 
even if ho and the redshifts zi, z s are given. 
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Table 1 

Plausible lens parameters to fit the images of a perfect Einstein cross and time delays with Ho = 72 km/s/Mpc 



Size 


Mass 


Shear 


Source Convergence 




R e / Re 




(71,72) 


(X S ,Y S ) k(0) 


k(R e ) k{10R e ) k(1Q0R e ) 


0.55/1.4 


1.6 


(-0.11,0.14) 


(-0.05,0.18) 5.3-1.6x10-^ 


1Q- 1 10~ 4 + 0.3^ 10~ 9 + 0.01i; 




Fig. 1. — The function N 4 (r) and the first and second derivatives N'(r) and N"(r). The images are at the rescaled radius r = 1. 
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Fig. 2. — Panel (a) shows time delay contours at intervals of 1 day for the four lens models. All models reproduce the same image positions 
and time delay with ho = 0.72. The difference of the models starts to be visible (the contours become thicker) beyond 3 arcsecs from the 
center. Panel (b) shows cuts of the time delay surfaces of the four lens models along the radial direction from the lens to the time delay 
minima (marked by "C" and "Al"on solid curves) and from the lens to the saddle images (marked by "B" and "A2" on dashed curves). From 
bottom to top lines are models with decreasing amount of halos v = 0.16,0.08,0.04,0 with color coded as red, blue, green and black. The 
measured delay tsc = 25 days implies ho = 0.72 for all these models. The time delay ratio tAcl^BA = 1-13 ±0.18 of Barkana (1997) is also 
reproduced. 
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R arcsec 



Fig. 3. — The deflection strength of the lens u 2 = = 44 as a function of the distance from the lens; this is effectively a rotation curve 

of the lens. From top right to down right, from thinner to thicker lines are models with decreasing amount of halos u = 0.16, 0.08, 0.04, with 
color coded as red, blue, green and black. Note that both Keplerian and flat rotation curve models are allowed by the data. 



